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I.  Summary  of  accomplishments 

1.  Adsorption  of  small  Au  andAg  clusters 

Experiments  have  shown  recently  that  small  Au  clusters  catalyze  a  number  of  reactions 
even  though  the  bulk  gold  is  chemically  inert.  We  have  shown1,  by  performing  a  number 
of  “computer  experiments”  that  the  activity  is  invariably  connected  to  surface  roughness. 
Recently  we  have  refined  this  picture2  to  show  that  it  is  the  geometry  of  the  LUMO  or 
HOMO  that  controls  activity,  not  the  geometry  of  the  cluster.  We  have  studied  adsorption 
of  Au  on  stoichiometric  Ti02  and  on  a  Ti02  surface  having  oxygen  vacancies.  We  have 
shown  that  the  electrons  left  behind  when  a  vacancy  is  formed  are  delocalized  around  the 
vacancy  and  not  localized  on  the  Ti  atoms.  Recent  STM  work  by  Makki  (ECOSS 
meeting,  Berlin,  2005)  confirms  this  picture.  We  have  also  shown  that  addition  of  Au  to 
the  vacancy  results  in  charge  transfer  to  the  Au  atom  and  this  causes  a  change  in  its 
chemical  activity.  It  is  very  strange  that  the  small  Aun  clusters  are  planar  and  we 

performed  high-level  calculations4  to  determine  the  value  of  n  for  which  a  three- 
dimensional  cluster  is  most  stable. 

2.  How  accurate  are  the  DFT  calculations? 

Density  functional  theory  is  the  only  quantum  chemical  method  used  to  perform 
calculations  on  systems  of  interest  to  catalysis.  We  need  to  know  what  kind  of  accuracy 
we  can  expect  from  such  calculations.  Since  it  is  very  difficult  to  characterize  precisely 
experimentally  catalytic  systems,  it  is  difficult  to  validate  DFT  by  comparison  with 
experiment.  We  have  chosen  to  evaluate  the  accuracy  of  DFT  by  applying  it  to  small 
clusters  for  which  we  can  perform  accurate  quantum  chemical  calculations.  We  found 
that  the  results  are  qualitatively  correct  for  oxygen  adsorption  on  small  Au  clusters5,  but 
they  are  accurate6  for  H2  and  propene  adsorption7.  In  the  latter  case  the  DFT  calculations 
also  agree  with  accurate  measurements  of  binding  energy  to  gas-phase  clusters, 
performed  in  Professor  Bowers’s  laboratory. 

3.  The  electronic  properties  of  linear  metal  clusters  deposited  on  a  metal  surface 
Wilson  Ho  has  performed  several  spectacular  experiments  whose  results  could  be 
interpreted  very  accurately  by  using  a  one-dimensional  model.  Unfortunately,  while  the 
data  was  fitted  well,  the  model  violated  Pauli  principle.  Our  theoretical  analysis8  showed 
that  the  one  must  use  a  three-dimensional  model  which  satisfied  Pauli  principle  and  also 
explained  why  the  experimental  result  behaved  as  if  the  system  was  one-dimensional. 

4.  Qualitative  rules  explaining  the  binding  of  molecule  to  Au  and  Ag  clusters 
We  have  developed  an  extremely  simple  model  which  showed  that  the  results  of  the 
density  functional  calculations,  for  binding  of  several  molecules  (CO,  02,  propene)  to  Au 
and  Ag  clusters,  can  be  predicted  by  examining  the  shape  and  the  energy  or  either  the 
LUMO  or  the  HOMO  orbitals9-11.  These  rules  have  been  used  recently  in  Bowers  group 
to  determine  the  structure  of  gas-phase  metal  clusters. 
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5.  A  proposal  of  a  new  system  for  catalytic  oxidation 

We  performed  density  functional  calculations  to  show  that  doped  oxides  might  provide  a 
new  class  of  oxidation  catalysts.  A  paper  on  this  subject  was  accepted  by  Catalysis 
Letters 12  and  two  others  are  being  prepared. 

6.  A  new  method  for  solving  scattering  and  bound  state  problems  in  quantum 
mechanics 

Plane  wave  basis  set  and  Fourier  transforms  are  extensively  used  in  quantum  mechanics 
for  solving  eigenvalue  and  dynamics  problems.  The  use  of  this  basis  set  has  several 
advantages:  accurate  representation  of  the  wave  function  and  rapid  and  very  accurate 
evaluation  of  derivative  and  Coulomb  integrals.  Unfortunately,  these  methods 
automatically  make  the  system  periodic  and  sometimes  this  causes  artifacts  in  the  physics 
of  the  system.  We  have  shown13-15  that  a  basis  set  of  Chebyshev  polynomials  has  all  the 
good  properties  of  the  plane  waves,  but  do  not  impose  periodicity.  We  have  also  shown 
that  this  basis  set  is  particularly  useful  when  one  intends  to  use  the  values  of  the  wave 
function  on  a  spatial  grid  as  the  unknown  quantities  in  the  calculation. 

7.  Is  phenomenological  kinetics  applicable  to  catalysis  problems? 

Phenomenological  kinetics,  of  the  kind  taught  in  the  introductory  physical  chemistry 
courses,  is  extensively  used  to  analyze  the  kinetics  of  the  catalytic  reactions.  The  results 
of  such  analysis  are  extensively  used  in  industry  for  reactor  design  and  in  science  for 
determining  the  reaction  mechanism.  Recently  the  kinetics  of  the  catalytic  reaction  has 
been  analyzed  by  using  the  kinetic  Monte  Carlo  method,  which  is  capable  to  treat 
realistically  all  the  microscopic  events  taking  place  during  a  reaction:  adsorption, 
migration  on  the  surface,  adsorbate  morphology,  reaction  and  product  desorption.  We 
have  performed16  extensive  comparisons  between  the  kinetic  Monte  Carlo  results  and  the 
phenomenological  theory,  for  the  reaction  of  CO  oxidation  by  Ru.  We  have  shown  that 
the  phenomenological  theory  is  qualitatively  correct  but  it  is  not  quantitatively  accurate. 
We  have  explained  and  analyzed  the  reasons  for  the  fai  lure  of  the  phenomenological 
approach. 
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